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Abstract. The Galactic synchrotron emission is expected to be the most relevant source of astrophysical contamination in 
cosmic microwave background polarization measurements, at least at frequencies y ^ 70 GHz and at angular scales ^ ^ 30'. We 
present a multifrequency analysis of the Leiden surveys, linear polarization surveys covering the Northern Celestial Hemisphere 
at five frequencies between 408 MHz and 1411 MHz. By implementing specific interpolation methods to deal with these 
irregularly sampled data, we produced maps of the polarized diff'use Galactic radio emission with a pixel size ^ 0.92°. We 
derived the angular power spectrum (APS) (PI, E, and B modes) of the synchrotron dominated radio emission as function of 
the multipole, £. We considered the whole covered region and some patches at diff'erent Galactic latitudes. By fitting the APS 
in terms of power laws (Q ~ k • i"), we found spectral indices that steepen with increasing frequency: from a ~ -(1-1.5) 

at 408 MHz to or (2-3) at 1411 MHz for 10 < ^ < 100 and from a ~ -0.7 to or ~ -1.5 for lower multipoles (the exact 

values depending on the considered sky region and polarization mode). The bulk of this flattening at lower frequencies can be 
interpreted in terms of Faraday depolarization eff'ects. We then considered the APS at various fixed multipoles and its frequency 
dependence. Using the APSs of the Leiden surveys at 820 MHz and 1411 MHz, we determined possible ranges for the rotation 
measure, RM, in the simple case of an interstellar medium slab model. Also taking into account the polarization degree at 
1.4 GHz, it is possible to break the degeneracy between the identified RM intervals. The most reasonable of them turned out 
to be RM ~ 9 - 17 rad/m^ although, given the uncertainty on the measured polarization degree, RM values in the interval 
~ 53 - 59 rad/m^ cannot be excluded. 
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1. Introduction 

In the last decade an impressive number of experiments has 
been dedicated to measurements of the cosmic microwave 
background (CMB) anisotropics, first observed in tempera- 
ture by CORE (Smo ot et al. 19901) . Modern cosmologies pre- 
dict also the existence of polarization fluctuations, produced 
at the recombination epoch via Thomson scattering (see, e.g., 
IKosowski 1999t) . However, the foreseen degree of polarization 
of the CMB anisotropics should be ^ 10%, implying a sig- 
nal much weaker than in temperature and therefore extremely 
difl&cult to reveal. The first detection of the CMB polariza- 
tion achieved by DASI ( Kovac 2002 ) and the recent measure 
by BOOMERanG JMasi et al. 20051 [Montroy et al. 2005| ) con- 
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firmed the expectations. The polarization data expected by the 
NASA WMAP ^ satellite will permit us to extend the informa- 
tion about polarization to a wider portion of the sky, though at 
small scales (^ ;^ 0.5°) its sensitivity will likely allow a detec- 
tion rather than an accurate measure. 

In the near future, the ESA Planck ^ satellite will signifi- 
cantly improve these results by observing the whole sky at nine 
diflTerent frequencies between 30 and 857 GHz, both in temper- 
ature and polarization, with unprecedented resolution and sen- 
sitivity (e.g. see ITaubeT 2004 and references therein). In par- 
ticular, the sensitivities per pixel of the Planck maps are ex- 
pected to be between 3-4 and 10 times better than in the WMAP 
data, producing a significant improvement in terms of angu- 
lar power spectrum recovery. Consequently, Planck alone will 
allow an estimate of the cosmological parameters much more 
precise than that obtained by WMAP and other CMB exper- 

^ http : //lambda. gsfc .nasa. gov/product/map/ 
^ http://www.rssd.esa.int/planck 
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iments combined with different cosmological data. Moreover, 
it will accurately measure the CMB £^-mode and could poten- 
tially observe the 5-mode as well (e.g. see |Burigana et al. 2004| 
and references therein). 

Strong efforts have been made to evaluate the impact 
of the foregrounds on the feasibility of CMB anisotropy 
measurements, as well as to understand how to separate 
the different components contributing to the observed 
maps. The methods elaborated for this purpose can be 
divided in two categories: blind (Independent Component 
Analysis, see Baccigalupi et al. 2000| IMaino et al. 2002 
Baccigalupi et al. 2004 Delabrouille et al. 2003 ) and non- 
blind (Wiener filtering and Maximum Entropy Methods, see 
respectively Bouchetet al. 19991 and iHobson et al. 1998). The 
latter class needs an a priori knowledge of the frequency and 
spatial dependence of the foreground properties, while a blind 
approach only requires maps at different frequencies as input. 
On the other hand, blind methods benefit from having ancillary 
data (such as realistic templates) to improve the separation 
robustness and quality ^. 

The Galactic polarized diffuse synchrotron radiation is ex- 
pected to play a major role at frequencies below 70 GHz 
on intermediate and large angular scales (^ ^ 30'), at least at 
medium and high Galactic latitudes where satellites have the 
clearest view of the CMB anisotropics. At about 1 GHz the syn- 
chrotron emission is the most important radiative mechanism 
out of the Galactic plane, while at low latitudes it is compa- 
rable with the bremsstrahlung; however the free-free emission 
is unpolarized, whereas the synchrotron radiation could reach 
a theoretical intrinsic degree of polarization of about 75%. 
Consequently, radio frequencies are the natural range for study- 
ing the Galactic synchrotron emission, although it might be af- 
fected by Faraday rotation and depolarization. Nowadays"^ the 
only available data suitable for studying the polarized Galactic 
synchrotron emission on large scales are the so-called Leiden 
surveys ( |Brouw & Spoelstra 1976| ). These are linear polariza- 
tion surveys, extending up to high Galactic latitudes, carried 
out at five frequencies between 408 MHz and 1411 MHz. 
We have elaborated specific interpolation methods to project 
these surveys into maps with a pixel size ^ 0.92°. These 
maps can be exploited in different contexts. They are suit- 
able to analyse the statistical properties of the Galactic ra- 
dio polarized emission. Beside this obvious application, they 
can be used as input to construct templates for simulation ac- 

^ Future component separation methods will probably search for a 
good compromise between the relaxation of the a-priori assumptions 
on the independent signals superimposed on the microwave sky and 
the exploitation of the physical correlation among the various compo- 
nents. Also it will be crucial to analyse microwave data together with 
ancillary data from radio and infrared frequencies. 

^ During the final phase of the revision process of 
this work, the DRAO 1.4 GHz polarization survey has 
been released ( Wolleben et al. 2006 1. The data are avail- 
able at http://www.mpifr-bonn.mpg.de/div/konti/26msurvey/ 
or http://www.drao-ofr.hia-iha.nrc-cnrc.gc.ca/26msurvey/. or 
^ffp ://cdsweb.u- strasbg.fr/cgi-bin/qcat/? J/A+ A/448/4 1 1^ See La 
Porta et al. 2006 for a first angular power spectrum analysis of the 
DRAO survey. 



tivities in the context of current and future microwave polar- 
ization anisotropy experiments. For example they can be uti- 
lized to build templates for straylight evaluation (see, e.g., 
IChallinor et ar2000.,Burigana et al. 200T1 IBarnes et al. 200311 
and for component separation analyses. 

It is standard practice to characterize the CMB anisotropics 
in terms of the angular power spectrum ^ (APS) as a function 
of the multipole, i (inversely proportional to the angular scale, 
i ^ 180/6>(°)). The APS is an estimator related to the angu- 
lar two-point correlation function of the fluctuations of a sky 
field JPeeb les 1993). The APS is not able to fully character- 
ize the complexity of the intrinsically non-Gaussian Galactic 
emission. Nevertheless, it is of increasing use also in the study 
of the correlation properties of the Galactic foregrounds (see, 
e.g., Baccigalu pi et al. 200 1| for an application to the polar- 
ized radio emission). We stick to the commonly adopted ap- 
proach and consider the APS of both the polarized intensity, 
PI, and the E and B modes ^ (see IKamionkowski et al. 19971 
and |Zaldarriaga 20"QT] ). 

In Sect. 2 we briefly summarize the properties of the Leiden 
surveys and of some selected areas relevant for the analysis at 
the smaller angular scales. In Sect. 3 we introduce the gen- 
eral criteria adopted to elaborate the class of our interpolation 
algorithms. The final, optimized version of the algorithm and 
the accuracy of the derived maps and APSs are described in 
Sect. 4. In Sect. 5 we present our results in terms of maps, 
APSs, and a parametric description of the APSs. A brief com- 
parison with a preliminary version of the DRAO polarization 
survey (Wolleben et al. 2003; Wolleben 2005 ) at 1.4 GHz is 
given in Sect. 6 in terms of APS. Sect. 7 is devoted to the mul- 
tifrequency analysis of our results and to their interpretation in 
the context of an interstellar medium slab model. We summa- 
rize and discuss our main results in Sect. 8. 

2. The Leiden surveys and selected areas 

The so-called Leiden surveys are the final results of different 
observational campaigns carried out in the sixties with the 25- 
m radio telescope at Dwingeloo. The complete data sets as well 
as the observation reduction and calibration methods have been 
presented by Brouw & Spoelstra (1976). The observations were 
performed at 408, 465, 610, 820 and 1411 MHz with an angular 
resolution of Ohpbw = 2.3°, 2.0°, 1.5°, 1.0°, 0.6°, respectively. 
At each frequency, the observed positions are inhomogenously 
spread over the Northern Galactic Hemisphere and globally 
cover a sky fraction of about ~ 50%. The errors, quoted for 
points observed at least twice, follow Gaussian distributions 
with mean values of 0.34, 0.33, 0.16, 0.1 1 and 0.06 K from the 
lower to the higher frequency, respectively. The observations at 
the same position were taken at different azimuths, elevations, 
and days (the adopted scanning strategy being random in coor- 
dinates and time). A robust control and removal of the ground 

^ It contains all the relevant statistical information in the case of 
pure Gaussian fluctuations. 

^ We keep here the physical dimension of the maps (antenna tem- 
perature, typically in K or mK) and consequently express the APS in 
K2 or mK2. 
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radiation contaminating the observations have been also per- 
formed. These properties imply a very low level of contami- 
nation from residual systematic effects, at least in comparison 
with the intrinsic survey sensitivities, and make these surveys a 
reference point (though insufficient by itself) for absolute cal- 
ibration of radio data. As an example, the absolute calibration 
of the Effelsberg maps at 1 .4 GHz has been performed using 
the Leiden survey at 141 1 MHz (see Uyanik er et al. 1 998 for a 
description of the procedure). 

The main problem in the analysis of the Leiden surveys is 
their poor sampling across the sky; it is necessary to project 
them onto maps with pixel size of about 3.7° (we use here the 
HEALPix scheme, see IGorski et al. 20051 to find at least one 
observation on each pixel of the observed region, that would 
limit the angular power spectrum recovery to ^ ^ 50. On the 
other hand, for some sky areas the sampling is significantly 
better than the average, by a factor ^ 4, and the maps appear 
filled even for a pixel size of ~ 1.8°, allowing one to reach 
multipoles ^ ^ 100. We identified three regions that permit one 
to investigate this multipole range at both low and middle/high 
Galactic latitudes: patch 1 [(110° < / < 160°, 0° < /? < 20°)]; 
patch 2 [(5° < / < 80°, b > 50°) together with (0° < / < 5°, 
b > 60°) and (335° < / < 360°, b > 60°)]; patch 3 [(10° < / < 
80°,/?>70°)]. 

We underline that in these patches the polarized signal 
and the corresponding signal-to-noise ratio are typically higher 
than the average; they are associated with the brightest struc- 
tures of the polarized radio sky, i.e. the "Fan Region" (patch 1) 
and the "North Polar Spur" (patch 2 and 3). The "North Polar 
Spur" (NPS) has been extensively studied (see lSalter 1983 and 
|Egger & Aschenbach 1995 J . The theories that better match the 
observations are all variants of the same model and interpret 
the NPS as the front shock of an evolved Supernova Remnant 
(SNR), whose distance should be ~ 100 + 20 pc (inferred 
from starlight polarization, see [Bingham 1967[ ). In contrast, 
the present knowledge of the "Fan Region" is much poorer. A 
rough estimate of its distance can be derived from purely geo- 
metrical considerations: it has such an extent on the sky that it 
must be located within (l-2)xl0^ pc from the Sun not to have 
an unrealistically large dimension. Wolleben et al. (2005) has 
recently proposed a diff'erent view of the "Fan Region": from 
the distances of HIT regions that they consider to be located in 
front of this structure and to depolarize its synchrotron emis- 
sion, they claim that the "Fan Region" is indeed an enormous 
object extending up to the Perseus Arm, i.e. up to kpc dis- 
tance. 

Another common and remarkable characteristic of the selected 
areas is that the polarization vectors appear mostly aligned at 
the two higher frequencies of the Leiden surveys. 

3. Map production algorithm and consistency tests 

A simple method to project the Leiden surveys into HEALPix 
maps is to average the observations falling in each pixel. 
However, the maps of the polarized intensity {PI) and of the 
Stokes parameters {Q,U) produced in this way show discon- 
tinuities on scales 6 ^ Opixd- These discontinuities tend to add 
spurious power in the recovered APS, particularly at multipoles 



£ ~ 180/^p/;c. In order to smooth these discontinuities we im- 
plemented a "class" of specific "interpolation" algorithms to 
generate maps with Opix ~ 0.92°, for the whole sky region 
covered by the surveys. To each pixel is assigned a weighted 
average of the signal values falling in its neighbourhood, i.e. 
YuiiPi ' Z Pi with Pi = l/(cr^ ' d^), where and ct/ are the 
signal and the error associated with the observation and di its 
angular distance from the pixel centre. The radius within which 
the average is computed varies pixel by pixel and is selected 
according to the following guidelines: (i) to have enough ob- 
servations (at least 3); (ii) to use only observations quite close 
to the considered pixel centre (less than a few degrees); (iii) to 
minimize the resulting average fractional change of the signal 
with the variation of the number of points and of the radius. 
The first two conditions imply that the interpolation is as local 
as possible and is performed using a reasonable number of data 
points. The third recipe is a convergency criterion. 

In order to check the reliability of the method, we gener- 
ated at each frequency several groups of PI, Q, and U maps, 
each corresponding to a diff'erent choice of the distance power 
in the average weight (l/(cr^ • J")) for various pixel sizes (rang- 
ing between 3.7° and 0.92°). A detailed description of the con- 
sistency tests together with preliminary results referring to a 
first version of the polarized intensity maps have been reported 
in La Porta (2001) and Burigana & La Porta (2002). In that 
work we produced simulated maps of white noise, expected to 
contamine the astrophysical signal contained in the maps, ac- 
cording to three diff'erent recipes. We compared the APS of the 
simulated noise and signal maps to identify the ranges of mul- 
tipoles statistically relevant (i.e. where the signal is not masked 
by the noise) for various considered sky coverages. We found 
that the analysis of the survey full coverage provides a good es- 
timate of the polarization APS for £ ~ [2, 50], while the better 
sampled patches (characterized by a higher S/N ratio on aver- 
age) allow us to investigate the multipole interval ~ [30, 100]. 
We emphasize that the current APS analysis has a statistical 
meaning and cannot provide details on the local properties of 
the synchrotron emission. 



4. Interpolation algorithm optimization and 
accuracy 

In order to evaluate and optimize the signal reconstruction 
quality provided by our algorithm, we performed a final test. 
Starting from simulated polarization maps having known sta- 
tistical properties, we created a table of data analogous to the 
Leiden survey ones and run our code to build PI, Q and U 
maps to be compared with the input skies. We focused on the 
1411 MHz case, because this is the most seriously undersam- 
pled both in terms of number of observations per beam and 
of observed points. Therefore, the maps produced at this fre- 
quency are potentially most influenced by the interpolation 
method efl'ects, particularly at multipoles larger than some tens. 
We assumed = = Q = /c • ^ ^, with values of the param- 
eters K and a in agreement with the mean law obtained for the 
Galactic polarized synchrotron emission in our first analysis of 
the 141 1 MHz survey (reported in |Burigana & La Porta 2002, ). 
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For e ~ [30,200] we adopted ^ k = \{K^) and a = -3. 
Applying the HEALPix facility Synfast to these angular power 
spectra, we generated maps of Q and U with a pixel size of 
^ 13.7' and a beamwidth Ohpbw - 36', as in the 1411 MHz 
survey, and used them to build the simulated table of data as 
described below. 

At each observed position of the Leiden survey 

a) We first assigned the Stokes parameters corresponding to 
the Synfast simulation pixel in which it would project and 

b) The corresponding rms given in the original table. 

c) We then added to each Stokes parameter a random value 
extracted from a Gaussian distribution in order to mimic 
the eff'ect of the instrumental white noise (see, e.g., 
[Burigana & Saez 2 003 ). The adopted standard deviation of 
this Gaussian distribution has been chosen equal to the 
mean error quoted for the considered original data. In this 
way the signal-to-noise ratio of the simulated table is sta- 
tistically very close to the real case one. 

We then obtained PI as ^/Q'^ -\- and calculated its error crpj 
according to the standard error propagation rules. 
Finally, we applied our interpolation method to the simulated 
data table and produced independently PI ,Q, and U simulated 
maps with Opix - 0.92°. We ran our code several times, chang- 
ing the weight n to be adopted in the average and the criteria 
for selecting the interpolation radius, so to sort out the best so- 
lution. 

As a first self-consistency test of the resulting maps, we veri- 
fied that typically the PI map can be considered equal to the 
^JQ'^ -\- within the uncertainty limit in 85% of the pixels. 
Then we compared the reconstructed maps with the input sky 
maps directly obtained with the Synfast facility (see Fig.[I}. 

We first degraded ^ the Synfast simulations to 0.92° pixel 
size maps and considered only the portion of the sky covered 
by the Leiden survey. Computing the APSs of both the original 
and the reprocessed maps we found that the best agreement 
between them is obtained: 

1. Adopting as polarized intensity the quantity 

-yjQreproc + ^reproc (inStCad of the PIreproc, produCCd 

applying the filling code to the PI values of the data table); 

2. Using the weight l/(cr^ • d^)', and 

3. Interpolating the data in the minimum circle containing at 
least two observations. 

The first choice is motivated by the necessity to (implicitly) 
preserve the information on the polarization angle and is a pos- 
teriori justified by the following result: the APS of PI map re- 
constructed applying the code directly on the table of PI values 
presents a certain loss in power with respect to the input case. 
The reason is likely that the PI interpolated value overestimates 
the eff'ective mean value of the field in the presence of changes 
in the signs of the Stokes parameters and consequently leads to 

^ For lower multipoles we assumed a slope < -2.5. Precisely, in 
this test we used a = -2.5 for £ ~ [6, 30] and or = for ^ < 6. 

^ We verified that this operation does not significantly aff'ect the 
APS on the multipole range of interest. 



an underestimation of the mean value of the field fluctuations. 
Criteria 2 and 3 imply that the interpolation is as local as pos- 
sible according to the sampling of each sky region. This ap- 
proach is quite diff'erent from that pursued by other authors 
dBruscoli et al. 20021 : they performed a Gaussian smoothing 
everywhere with constant Ohpbw (of 3°), neglecting the sam- 
pling dependence on the position in the sky and on the fre- 
quency. 

In principle, before computing any average we should have 
transformed the polarization vectors using the formula of the 
parallel transport on the sphere; however operating only a lo- 
cal interpolation, the consequent mixing of the two orthogonal 
components of the polarized signal turned out to be negligi- 
ble ^ 

Fig. |2| shows that the statistical information contained in the 
synthetic data table is appreciably traced out in the multipole 
ranges of interest both for the half- sky map and for the patches 
(for conciseness we report here only the results found for patch 
2). 

Fitting the APS of the reprocessed map, we quantified the pre- 
cision level at which the input values of the amplitude and of 
the spectral index were recovered: by considering several real- 
izations for the various cases, we found relative errors within 
^ 10% for both the parameters. 

We also checked the reliability of the maps themselves in a 
more quantitative way. The distributions of the pixel to pixel 
diff'erence between the two series of maps (i.e. PI, Q and U pro- 
duced by Synfast and reprocessed by our interpolation code) 
are, in all cases, close to Gaussian distributions, with standard 
deviations comparable with the typical mean error of the data 
table (see Table [I] and Fig.|3li. 

We repeated the test described above skipping the step c) 
The distributions of the diff'erences between the input and re- 
processed maps are again close to Gaussian distributions, but 
with standard deviations slightly smaller than those found in- 
cluding step c) (in particular, for the full survey coverage the 
standard deviation is now smaller by about 20%). This indi- 
cates the existence of a geometric contribution in the error dis- 
tribution over the maps related to the poor sampling of the data. 



^ In order to drag vectors on the sphere preserving strength and an- 
gles, one should apply the parallel transport formula estabilished by 
diff'erential geometry. Integrating on paths close to geodetic, a suit- 
able approximation is given by 
pf^^ = Pe ■ cos(0 + ■ sin(0 
P^cp^^ = -Pe ' sin(0 + Pep ■ cos(0 

where ^ = cos(6ave)'^(p Oave = (0in+0fin)/2 f scc lBruscoli et al. 20021 
for further details). Using these equations it is possible to estimate the 
percentage contribution of each initial Stokes parameter to the new 
parameter expression, once it has been moved to another point of the 
sphere. We assumed S(p = 3°, corresponding to the mean value of 
the interpolation radius in the selected patches. Neglecting the parallel 
transport in our interpolation method causes a mixing of the Stokes 
parameters of ^ 1% in the equatorial area and of ^ 5% in the medium 
and high latitude regions. 

Formally, in this case we apply the step c) but using in the noise 
generation a standard deviation negligible with respect to the values 
of P/, Q and U (the implemented code assumes that a mean error is 
given for each measurement of the observed quantities). 
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Fig. 1. Synfast realizations (left side) at 1411 MHz versus maps reconstructed applying the final version of the filling code to the 
synthetic data table (right side). See text for details. 



Expressing the rms of the pixel to pixel diff'erence distribution 
as 



map 



^samp ^WN 



the analysis of reprocessed maps (with and without noise) gives 

(TwN - const X (Tsamp 

with const ^ 0.63,0.87,0.90, 1.15, respectively for the whole 
survey coverage and for the patches 1, 2, and 3, without signif- 
icant diff'erences among the maps of PI, Q or U. 
The relative weight of the geometrical contribution to the over- 
all error decreases with the goodness of the sampling. At the 
other frequencies of the Leiden surveys the sampling and the 
signal-to-noise ratio are better than at 1411 MHz; as a conse- 
quence, we can look at the ratio between (Tsamp (resp. (Tmap) 
and the signal at 1411 MHz as a conservative upper limit to 
the relative error caused by geometric eff'ects (resp. geometric 
eff'ects plus experiment white noise) at the lower frequencies. 



To verify whether the goodness of the signal reconstruction 
depends on the input sky characteristics, we repeated the 
entire test adopting a diff'erent input APS, more appropriate at 
408 MHz. At this frequency the analysis of the preliminary 
polarization maps led to diff'erent results for the APSs than 
the 1411 MHz case. In particular, the spectral index was 

flatter (roughly 1.5 in the multipole range [30,200]), both 

for the whole coverage map and for the patches. The same 
test repeated with this slope confirmed the results described 
before. The quality of the APS reconstruction is slightly less 
accurate than in the previous case, however the relative errors 
estimated for the recovered amplitudes and spectral indices 
remain within ^ 10%. 




Multipole (/) Multipole (/) 



Fig. 2. Comparison between the APSs of the input maps (dot Hnes) extracted from the Synfast simulations (degraded to 0.92° 
pixel size maps) and the reprocessed maps (solid lines). From the bottom: • 10 and • 100. Note the excellent agreement 

between the input and reprocessed APSs up to ~ 100. 
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Table 1. The pixel to pixel difference between the maps directly generated with Synfast (after degradation to 0.92° pixel size 
maps) and the reconstructed maps contaminated by white noise, for the whole survey and the considered patches. Columns 2 and 
3 give the average and rms of the Gaussian function that best fits the diff'erence distribution (shown in Fig.O. Columns 4 and 5 
contain the average and standard deviation of the diff'erence distribution. Columns 6, 7 and 8 give the percentages of pixels with 
absolute diff'erence within the reported levels (here Icr corresponds to the rms reported in column 5). 



5. Results 



In Fig.|4|we show the polarization maps (PI, Q, and U) pro- 
duced at all frequencies with the final version of the code de- 
scribed in the previous section. 

Using the Anafast facility of the HEALPix package, we 
computed the APSs of the polarized intensity and of the E and 



B modes for the half-sky coverage and for the patches at the 
five frequencies of the surveys. In addition to the three better 
sampled areas we considered another region (patch 4) located 
at (70° < / < 120°, -45° < b < -15°). This patch has a sky 
sampling in the average (its APS is then statistically relevant 
only for ^ ;^ 50) and is characterized by a rather low signal. 



The monopole of each map has been subtracted before APS esti- 
mation and the obtained Q have been divided by the fractional cover- 
age of the considered area to renormalize it to the whole sky case. 
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Fig. 3. Distribution of the pixel to pixel differences between the maps directly simulated with Synfast (degraded to 0.92° pixel 
size maps) and the reprocessed maps at 141 1 MHz. We display here the results of the test, both with (upper panels) and without 
(lower panels) white noise contamination, in the case of the full coverage 1411 MHz Leiden survey. The superimposed curves 
are the Gaussian distributions giving the best fit to the data. 
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Table 2. Values of the best fit parameters k and a obtained in the analysis of the survey full coverage APS at all the frequencies 
(as shown in Fig|6|). Relative errors on the best fit parameters are ~ 10%. Two multipole ranges are considered: [2, 1] and [^, 200] 
was also a free parameter in the fit). EB refers to (C^ + C^)/2. 
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Table 3. Values of best fit parameters k and a obtained in the analysis of the APS of the patches at all the frequencies (as shown in 
FigQ. Relative errors on the best fit parameters are ~ 10%. The considered multipole range is [30-200]. refers to (Cf -\-C^)/2. 



The E and B modes turned out to be extremely similar, so their average (C^ + C^)/2. In some cases (namely for patch 1 
that, instead of studying each of them singularly, we consider and 2 at 1411 MHz and patch 2 and 4 at 408 MHz) we could 
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Fig. 4. PI, 2, and U maps obtained from the Leiden surveys data tables using the final version of the interpolation algorithm at 
all frequencies. 



even consider a unique APS defined as (C^^ + + C^)/3. 
The maps represent the Galactic polarized synchrotron emis- 
sion, smoothed with the beam of the radiotelescope and con- 
taminated by the noise. As a consequence, their APSs can be 
fit as sum of two components, Q = C^^^^^Wt -h C^. We ex- 



ploit the power law approximation C 



synch 



and as- 



sume a symmetric, Gaussian beam, i.e. a window function 
= where (Tb = ^//pw(rad)/ V81n2 . Under the 

hypothesis of uncorrected Gaussian random noise (white 
noise), = ~ const . We have performed the fit on the 



Even in the ideal case of negligible systematic effects in the data, 
the averaging and the interpolation over the poorly sampled data may 
imply deviations from this simple assumption. Moreover, the noise 
level in the maps will be connected not only to the instrumental noise, 
but also to the sampling. A posteriori, the APS flattening found at high 
£ justifies our simple approximation. 



multipole range [30, 200] for the patches and [2, 200] for the 
full-coverage maps. In fact, the flattening occurring at higher 
multipoles helps recover the noise constant, in spite of the fact 
that no reliable astrophysical information can be derived for 
^ > 100 from the APSs even in the three better sampled re- 
gions ( |Burigana & La Porta 2002 |l. For the patches, we found 
that a single set of the parameters /c, a allows us to describe the 
APS in the entire range of multipoles. In the case of the survey 
full coverage we needed two sets of k, a, each of them appro- 
priate to a certain multipole range. At all frequencies a change 
in the APS slope occurs corresponding to ^ ~ 10 for both C^^ 
and (Cf Cf )/2. 

As an example, we display in Fig. 13 the APS and the corre- 
sponding best fit curve at 610 MHz for the survey full coverage 
and one patch. We show in Fig. |6|(resp. FigQ the recovered 
synchrotron term, Cf^^^, and list in Table |2|(resp. Table [Sj the 
corresponding k, a obtained for the survey full coverage (resp. 
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Fig. 5. Polarization APSs at 610 MHz (solid lines) for the survey full coverage (left panel) and patch 3 (right panel) together 
w^ith the corresponding best fit curves (dot-dash lines). The lower curves in each panel are C^^, while the other ones represent 
(Cf + Cf )/2 multiplied by 10. 

Polarized intensity E and B modes 
10^ F ' ' ' ' ' ' ^ ' ' ' ' ' 




Fig. 6. Synchrotron component of the APS derived for the survey full coverage for the polarized intensity (left panel) and for 
(Cf + C^)/2 (right panel) at the various frequencies: 408 (solid line), 465 (dashes), 610 (dot-dash), 820 (three dot-dash), and 
1411 MHz (long dashes). 



for the patches). In the multipole range 30 ^ £ ^ 100 we found 
a general trend: from the lowest to the highest frequency the 
slope increases from (1-1. 5) to (2-3), with a weak de- 
pendence on the considered sky region. 

6. Comparison with the DRAO polarization survey 

Soon, a new all-sky survey at 1.4 GHz will become available. It 
results from the merging of two diff'erent polarization surveys 
having the same angular resolution of about 36': the Northern 
part of the sky has been observed with the DRAO 26 m tele- 
scope ( Wolleben et al. 2003 ) . for the Southern hemisphere the 
30 m Villa Elisa telescope was used dTestori et al. 2003ft . The 
first was meant to compensate for the lack of information typ- 
ical of the Leiden surveys and it improves their sampling and 
sensitivity. However the Leiden surveys remain so far a unique 
tool for absolute calibration in the Northern sky, therefore the 
DRAO data have been tightened to them. A second observing 



period at DRAO has recently been completed. The final map is 
fully sampled in right ascension and has a spacing in declina- 
tion ranging between 0.25° and 2.5° (IWolleben et al. 200611 . In 
Fig. [HI we compare the APS of C^^, Cf, and derived from 
a preliminary version of the DRAO survey with the ones of the 
maps we constructed using the Leiden survey at 1411 MHz. 
For the full coverage maps the APSs almost coincide at ^ ;^ 10. 
At ^ ^ few X 10 the slopes are in good agreement in all cases, 
while the amplitudes derived from the Leiden survey are larger 
by a factor ~ 2.5 for the full-coverage maps and ~ 4 for the 
patches. As a consequence, the fitted amplitudes reported in the 
previous section might translate in an overestimate of the syn- 
chrotron emission APS at 1.4 GHz. We checked that the signal 
in the DRAO survey is typically lower than in map we con- 
structed from the Leiden survey at 1411 MHz, which seems 
consistent with the amplitude discrepancy found. A detailed 
analysis of the final version of the recent 1.4 GHz polariza- 
tion data and its comparison with the Leiden 1 .4 GHz survey 
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Fig. 7. As in Fig.|6lbut for the selected patches: patch 1 (top panels), patch 2 (middle panels), and patch 3 (bottom panels). 



is in progress. We stress here the good agreement between the 
slopes of the APSs obtained from the two data sets that con- 
firms the steepness of the APS we found at 1 .4 GHz from the 
Leiden survey. 

7. Multifrequency analysis 

In Sect. 5 we studied the APS of the Galactic diff'use polarized 
emission at each frequency, focusing on its multipole depen- 
dence. We first consider here the behaviour of the synchrotron 
APS at some fixed values of the multipole as function of fre- 
quency, being aware of the depolarization phenomena. In the 
second subsection we provide a qualitative interpretation of the 
APS steepening with increasing frequency in terms of Faraday 
depolarization eff'ects. 

7.1. Frequency dependence of the APS amplitude and 
Faraday depolarization 

The theoretical intrinsic behaviour of the synchrotron total in- 
tensity emission predicts a power law dependence of the APS 
amplitude on frequency: 

The values of JS deduced from radio observations change with 
the considered sky position and range between ~ 0.5 and 
~ 1 (see Reich & Reich 1988 and Platania et al. 1998). 
In principle, one would expect the synchrotron emission to 
be highly polarized, up to a maximum of ~ 75% (see 
[Gmzburg & Syrovatskii 1965| l. In the ideal case of constant de- 
gree of polarization the above formulae would apply also to 
the polarized component of the synchrotron emission, once the 
brightness temperature has been properly rescaled. However, 
depolarization eff'ects are relevant at radio frequencies, so that 
the predicted correlation between total and polarized intensity 
do not show up at all (except for the NPS and the Fan Region 



that are clearly evident both in total intensity and polarization 
maps at about 1.4 GHz). The current knowledge of for the 
Galactic diff'use polarized emission is very poor, but already 
indicates that due to depolarization phenomena the observed /3 
can be much lower than ~ 0.7 (e.g. for the NCP Vinyajkin & 
Razin 2002 quote /3 ~ -0.13). 

In Fig. 121 we plot C^^^^^(v) at boundary and intermediate 
values of the statistically significant multipole range in each 
coverage case. For the patches we chose £ = 30,50,70, 100, 
whereas for the half-sky maps we considered ^ = 5, 30, 50, 70. 
The APS amplitude actually decreases with frequency, but our 
results are not well fitted by a power law and exhibit shapes sig- 
nificantly flatter than that of the observed synchrotron emission 
in total intensity (approximately Cf^^^{v) oc y"^ "^), as expected 
in the presence of frequency-dependent depolarization efl'ects. 

An electromagnetic wave travelling through a magnetized 
plasma will undergo a change in the polarization status. In par- 
ticular its polarization vector will be rotated by an angle 

A0[rad] = 7?M[rad/m^] • A^[m^] . 

Here the rotation measure, RM, is the line of sight integral 

where Ue is the electron density and B\\ is the component of 
magnetic field along the line of sight. 

Given A0 oc y"^, the lower the frequency the greater the 
importance of the Faraday rotation eff'ect. Approximating the 
interstellar medium (ISM) as a superposition of infinitesimal 
layers each characterized by constant rie and B\\ (slab model, 
see Burn 1966), the observed polarization intensity relates to 
the intrinsic one according to the expression: 

TP^^^'(A) = TP^^^'^A) • |sinA0/A0| . 

Therefore the net eff'ect of a diff'erential rotation of the polariza- 
tion angle is a depletion of the polarized intensity, commonly 
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Fig. 8. Comparison between the DRAO APSs (solid lines) and the APSs of the maps reconstructed by our code (dashes) from 
the 1411 MHz Leiden survey properly rescaled according to Wolleben (2005). In each panel, the three pairs of curves (from the 
bottom to the top) correspond to the PI, E, and B modes; the E and B APSs have been multiplied respectively by 20 and 700. 



referred to as Faraday depolarization. Given (r^)^ oc Q(y) the 
above formula reads: 

Q(y) = Cf» • (sin A0/A0)^ 

where 0^^^(y) = k - £^ - y^^"^"^^ . At a given multipole, the ratio 
between the APSs at two frequencies can be written as: 

Q(yi)/Q(y2) = (yi/y2)-'^-'-(A02/A0i)'-(sinA0i/sinA02)' • 

Given a certain value of fi, the expression on the right side 
of this equation can be displayed as a function of RM (see 
Fig. [Toll . The RM intervals in which the curve (obtained for 
a certain JS) intersects the constant line representing the ob- 
served value of this ratio provide information about the pos- 
sible RM values This holds the in case of pure Faraday de- 
polarization. In reality there are at least two other phenomena 
affecting the observed polarized intensity, the beamwidth and 
the bandwidth depolarization (see IGardner & Whiteoak 1966 
and Sokoloff et al. 1998 ). If the polarization angles vary within 
the beam area, then observations of the same sky region with 
different angular resolution will give an apparent change of the 
polarized intensity with frequency (beamwidth depolarization). 
Suppose 9hpbw,i and 6hpbw,2(< Ohpbw,i) are the angular reso- 
lutions respectively at vi and y2. The effect can be removed 
by smoothing the map at the frequency y2 to the (lower) angu- 
lar resolution of the map at the frequency vi or, equivalently, 
multiplying the observed APSs ratio, Q(yi)/Q(y2) by e^^^^^^ , 
where Acr^ = cri^ - 0-2^ ((Ti = Ohpbwj (rad)/ V81n2 ; / = 1,2). 
This correction has been applied to convert the APSs shown 
in Fig.|9|to the ratios shown in Fig.fTOl 

Clearly, the RM values identified with this method do not refer to 
a precise direction in the sky, but to the considered sky area. 

Once the 1411 MHz map is properly smoothed to the lower res- 
olution of the 820 MHz survey, the bulk of possible differencies due 



Another depolarization effect is introduced by the receiver 
bandwidth: the polarization angle will change by an amount 

Al// = -RM • Ao^ • 2Ay/yo , 

where Aq and vq are the values at the bandwidth centre. The po- 
larization of the incoming radiation will be reduced by a factor 
sin Ai///Ai//, that in our case is negligible (Ay =1.7 MHz at the 
two lower frequencies. Ay = 3.6 MHz at the intermediate ones 
and Ay = 1.2 MHz at the highest one). 

In Fig.Iiniwe plot the C'f'\yi)/C'f'\y2) dXl = 30,70 
(resp. at ^ = 5, 30) in the case of the patches (resp. in the case 
of the full-coverage maps). 

The information coming from the lower frequencies is too 
poor to be able to learn something about 7?Ms; on the other 
hand, the ratio between the APSs at 820 and 1411 MHz give 
some indications. For a reasonable range of synchrotron emis- 
sion frequency spectral indices, the following values are com- 
patible with the data (see left-bottom panel of Fig.[I2l>: the rel- 
atively spread intervals 9-17, 53-60, 75-87, 123-130 rad/m^ to- 
gether with two narrow intervals at 70 and 140 rad/m^. In gen- 
eral there is a rather good agreement among the results coming 
from the patches and the full-coverage map. Spoelstra (1984) 
reported as a typical value RM = 8 rad/m^, however this is 
likely a lower limit to the real values (in addition, beam de- 
polarization effects were not taken into account in that study, 
likely implying a certain underestimation of RM). 

Another constraint to the possible values of the RMs comes 
from the polarization degree, defined as 

H = TP/T . 



to beam depolarization effects is eliminated and no longer contributes 
to the value of the ratio between the APSs at the two frequencies. 
Consequently, beam depolarization is not an issue in the following 
analysis. 
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Fig. 9. Synchrotron PI APS at all frequencies for some fixed representative multipoles, I close to the boundary and in the middle 
of the range of validity of our analysis. We considered £ = 30, 50, 70, 100 for the patches and ^ = 5, 30, 50, 70 for the full coverage 
maps. The diamonds represent the C^J^^^{v) at the survey frequencies. Diff'erent lines correspond to diff'erent sky regions, namely: 
full survey coverages (solid line), patch 1 (dots), patch 2 (dashes), patch 3 (dot-dashes), patch 4 (three dots-dashes). 



Exploiting the total intensity northern sky map at 1.4 GHz 
(Reich & Reich 1988), we estimate the mean polarization de- 
gree n in the portion of the sky covered by the Leiden survey 
at 1411 MHz. We find ff^' ~ 17, 27, 28, 25, 12% for the full 
coverage and for the patch 1, 2, 3, and 4, respectively. By us- 
ing the preliminary DRAG survey, for the same sky regions 

we obtaine slightly lower values, i.e. ^ ~ 12, 22, 24, 20, and 
9%, respectively. Given the relation existing between the intrin- 
sic and observed brightness polarized temperature in the slab 
model approximation, the observed degree of polarization can 
be written as: 

n^^^=ff"'^-|sinA0/A0|. 

In Fig.[n]the maximum possible value of H^^^ (derived assum- 
ing n^"^^ = 0.75 ) is displayed as a function of the RM. Each 
horizontal line represents the mean polarization degree actu- 
ally observed for the Leiden surveys in the considered coverage 
cases^^. 

Note that the observed mean polarization degree might be due 
to the combination of Faraday depolarization and beam depolariza- 
tion. If we were able to correct for beam depolarization effects, these 
lines would be shifted upwards in Fig.[n] leading to smaller values of 
RMs. Therefore the results of our analysis are rather conservative and 
provide upper limits to the RMs compatible with observations in the 
considered patches. 



Except for patch 4, the observed polarization degree sets 
an upper limit to the RM of ~ 50 rad/m^ (having assumed in 
Fig. [H] the maximum degree of polarization). Consequently 
the lowest of the RM intervals previously identified (exploit- 
ing the APSs amplitude ratio at the two highest frequencies of 
the Leiden surveys) seems to be the most reasonable. 

The above considerations hold under the following assump- 
tions: /) the physical properties of the magneto-ionic ISM re- 
sponsible for depolarization are homogeneous along the line 
of sight; //) the bulk of the observed polarization signal comes 
spatially from the same regions at all frequencies. Under sig- 
nificant violations of these assumptions the meaning of derived 
RM is by definition unclear. 

The real situation is generally much more complicated, due 
to the existence of turbulence in both the Galactic mag- 
netic field and the ISM electron density depolarizing the 
diff'use synchrotron emission (see Haverkorn et al. 2004). 
Another aspect is the presence of foreground magneto-ionic 
structures (Faraday screens) modulating the background syn- 
chrotron diff'use emission (see IWolleben & Reich 20041 and 
Uyaniker 2003 ) , which renders the understanding of the lo- 
cal ISM very difficult. There is a lack of information about the 
nature and spatial distribution of these structures and, conse- 
quently, about the corresponding depolarization eff'ects. 
Nevertheless, in patch 1, 2 and 3 the polarization vectors are ex- 
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Fig. 10. Ratios (C~^"^ (yi)/C~^"^ (y2)) between the polarized synchrotron APSs at two frequencies for a pair of representative 
multipoles in the range of validity of our analysis. For the patches we plot the ratio at ^ = 30 and 70, whereas for the full coverage 
maps we chose 1 = 5 and 30. The horizontal lines represent the observed ratio corrected for beamwidth depolarization; thick 
lines refer to the lower of the two considered multipoles. Different sky regions are considered: full survey coverages (solid line), 
patch 1 (dots), patch 2 (dashes), patch 3 (dot-dashes), patch 4 (three dots-dashes). The oscillating curves represent the theoretical 
ratio as a function of RM in case of pure Faraday depolarization for three choices of the synchrotron emission frequency spectral 
index: -(p + 2) = -2.5 (dashes), soHd -(]3 + 2) = -2.75 (soHd lines), and -(]3 + 2) = -3 (dot-dashes) . The last panel show the 
same three curves in the hypothetical case of availability of observations at higher frequencies. 



0.8 




Fig. 11. The maximum possible value of the observed polarization degree (i.e. = IF^^^ • | sin A0/A0I, where IF^^^ = 0.75) 
is displayed as function of the RM. The horizontal lines represent the mean polarization degree observed in the considered sky 
regions. Survey full coverage (solid line), patch 1 (dots), patch 2 (dashes), patch 3 (dot-dash), patch 4 (three dots-dash). 
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tremely well aligned and similarly distributed at the two higher 
frequencies of the Leiden surveys. This observational fact sug- 
gests that for these areas the slab model is a reasonable approx- 
imation of the local ISM, at least for the angular scales inves- 
tigated in the present analysis. Consequently, the above con- 
siderations on the 7?Ms derived from the the APS ratios should 
be basically correct with respect to these regions (whereas they 
should be taken with caution in the other considered cases). 

Finally, the right-bottom panel of Fig. fTHI shows, as an ex- 
ample, that a multifrequency sky mapping at v ~ some GHz 
would provide direct information on Faraday depolarization, 
significantly reducing the degeneracy on RM (unavoidable in 
maps limited to v ~ 1 GHz). 

7.2. Effect of the Faraday depolarization on the APS 
slope 

Our analysis of the Leiden surveys showed that the synchrotron 
emission APS steepens with increasing frequency (see Sect. 5). 
One possible explanation for such a behaviour of the APS re- 
lies on depolarization arguments: fluctuations of the electron 
density and/or of the magnetic field (in strength and/or direc- 
tion) might redistribute the synchrotron emission power from 
the large to the small angular scales, creating fake structures. 
As a consequence the corresponding APS will become flat- 
ter, the efl'ect being more important at the lower frequencies. 
Sticking to the slab model to describe the ISM, we simulated 
the Faraday depolarization efl'ects on the APS at 408 MHz 
and 1420 MHz as foUows. We assumed that the 1420 MHz 
polarization map JWoUeben et al. 2006t represents the intrin- 
sic Galactic synchrotron emission and used a spectral index 
of -2.75 to scale it down to 408 MHz^^. The polarization an- 
gle map at the two frequencies should be the same in the ab- 
sence of depolarization efl'ects; we adopted the polarization an- 
gle map at 1420 MHz as intrinsic. These four maps are the 
input skies for our toy- simulation of Faraday depolarization. 
To build an RM map we exploited the catalog produced by 
Spoelstra (1984)^^ ; it contains ~ 1000 RM data, whose min- 
imum and maximum angular distance is dmin ~ 0.02° and 
dmax ~ 9°. We interpolated these data to generate an RM map 



The spectral index we used is rather steep for the polarized emis- 
sion. However the present exercise is meant to justify the change of the 
APS slope with increasing frequency. It is not intended to reproduce 
the 408 MHz map reconstructed from the Leiden data. This would 
require at least a frequency spectral index map appropriate to the po- 
larization signal. The current knowledge of the spectral index distri- 
bution is limited to the total intensity radio diffuse Galactic emission 
and cannot be applied to the polarized emission, since the two come 
from different spatial regions. 

The RM map obtained interpolating the values of Spoelstra (1984) 
are most likely not reliable. Those values were obtained performing 
a linear fit to the five frequencies data of the Leiden surveys and 
assuming as a result of the fit the lower possible RM. Furthermore 
beam widths effects were considered to be negligible. Consequently 
the RMs derived by Spoelstra (1984) should be taken as indicative 
values. We simply exploited the existence of these data to perform toy- 
simulation of differential Faraday depolarization effects on the APSs 
of polarized emission. 



with Opixei - 1.84° (see Fig.FT^: the algorithm searches in 
circles of increasing radius from each pixel centre, considering 
radii from dmin to dmax- It stops as soon as it finds at least 3 
points and associates to the pixel the weighted average of the 
corresponding values; the weight is inversely proportional to 
d'^, where d is the distance from the pixel center. The results 
we show below have been obtained choosing n = I, although 
we checked that they remain valid also for n = 0.5, 2, 4. 
We use the RM map to transform the intrinsic and PI for 
Faraday rotation and differential Faraday depolarization: 

(poutput ~^ (pinput + ^0 
Ploutput PUnput • |^m(A0)/A0| 

. At 1420 MHz the input and output maps are very similar. 
At 408 MHz the situation is completely different: the large- 
scale structure appears significantly depressed in the depolar- 
ized map, whereas the small scale features become much more 
relevant (see Fig.fT2l. 

In Fig. [21 we compare the APSs of the input and the out- 
put maps at both the frequencies. At 1420 MHz the APSs are 
almost identical; on the contrary at 408 MHz the depolarized 
map APSs significantly flatten, as expected from the increase 
of small scale features. This supports the interpretation of the 
bulk of the APS flattening with decreasing frequency resulting 
from our analysis of the Leiden surveys in terms of the Faraday 
depolarization effect. 

8. Discussion and conclusions 

The Leiden surveys have been up to now the only available 
data at y ~ 1 GHz suitable for a multifrequency study of the 
polarized diffuse component of the Galactic synchrotron emis- 
sion on large angular scales. A linear polarization survey of the 
Southern sky had previously been carried out at 408 MHz us- 
ing the Parkes telescope (Mathewson & Milne 1965 ); unfortu- 
nately these data have never been put in digital form. In Fig. [21 
an image of the data from this survey has been combined with 
the Leiden survey at 408 MHz to show an all- sky linear polar- 
ization map. 

The "North Polar Spur", the "Fan Region", and the "Cetus 
Arc" (see Berkhuijs en 1971| for a sketch of the Galactic Loops 
these structures are thought to belong to) are visible as areas 
where the polarization vectors are longer (the length is propor- 
tional to the polarized intensity) and more regularly aligned, 
drawing small circles in the sky. The most plausible theory sug- 
gests that these Galactic spurs are the brighter ridges of evolved 
Supernova Remnant sitting relatively close to the Sun (within a 
few xlO^ pc) and obscuring our view of the Galactic diffuse 
non-thermal emission. The APSs obtained for these regions 
will be characterized by an amplitude considerably higher than 
the average, representing therefore an upper limit for the polar- 
ized signal of the Galactic synchrotron radiation. On the other 
hand, the APS should be steeper in these regions, since the 

The number of data points is about half the lowest one of the 
Leiden surveys catalogs, therefore we could not work with Opi^ei - 
0.92°. 
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Fig. 12. Maps of the Stokes parameters before (upper panels) and after (lower panels) the transformation for Faraday depolariza- 
tion. In the middle is placed the RM map used in the simulation. All maps have been smoothed to 2.3°. In the depolarized maps 
the large scale structure appear significantly depressed respect to the input ones, whereas the small scale structures are more 
relevant. 
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Fig. 13. E and B modes APSs before (solid lines) and after (dashed lines) the transformation according to Faraday depolarization. 
The upper panels refer to 1420 MHz and the lower ones to 408 MHz. The diff'erent importance of the depolarization eff'ects at 
the two frequencies is evident: the 1420 MHz APSs are almost unaff'ected, whereas the 408 MHz APSs significantly flatten. 
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Fig. 14. Map of polarization vectors at 408 MHz, obtained combining the corresponding figures of the Leiden survey (Brouw 
& Spoelstra 1976) and of the Parkes southern sky survey dMathev^son & Milne 1965 ). Unfortunately, the vector length has a 
diff'erent scale in the two original figures, therefore only a direct comparison of the direction is feasible. 
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polarization vectors appear rather ordered and with approxi- 
mately constant length, implying a relatively low fluctuation 
power at small scales. These considerations should be kept in 
mind when extrapolating the results to the microwave range in 
order to check the impact of the Galactic synchrotron emission 
on CMB polarization measurements. 

Implementing an interpolation method to deal with the sparse 
and irregular sampling of the Leiden surveys, we produced 
maps of PI, Q and U with 6 pixel - 0.92° at five frequencies 
between 408 and 1411 MHz. We accomplished several tests in 
order to check the reliability of the interpolation algorithm and 
to optimize its quality (see Sects. 3 and 4). We emphasize that 
the APSs computed from our maps are significantly diff'erent 
from those obtained by other authors (e.g. lBruscoli et al. "2 002) 
for similar areas of the sky, using the same data; namely, at 
the higher frequencies of the Leiden surveys we find steeper 
APSs (first recognized in Burigan a & La Porta 2002| l. A firm 
confirmation of our approach comes from the comparison of 
the APSs obtained for the 1411 MHz map with those of a pre- 
liminary version of the DRAO 1420 MHz survey, characterized 
by much denser sampling and higher sensitivity. 

We used the interpolated full coverage maps to study the 
APS of the polarized synchrotron emission on the multipole 
range [2,50], being limited at high £ by the map sensitivity 
(related to the accuracy and sampling of the original data). 
Exploiting some better sampled regions, approximately coin- 
ciding with the two most prominent Galactic spurs, we could 
extend our analysis up to ^ ~ 100. 

At each frequency the APSs of the entire coverage maps 
turn out to change slope at ^ ~ 10, C^^, C^, and presenting 
a steepening going from the lower to the higher multipoles (see 
Sect. 5). A general and remarkable result (see Sect. 5) is that 
the slope of the polarized synchrotron emission APS increases 

from 408 MHz to 1411 MHz; it changes from a (1-1.5) to 

a ~ -(2-3) for ^ ~ 10 to ~ 100 and from a ~ -0.7 to a ~ -1.5 
for lower multipoles, the exact value depending on the consid- 
ered sky region and polarization mode. 

One possible explanation for such a behaviour of the APS re- 
lies on Faraday depolarization arguments. Fluctuations of the 
electron density and/or of the magnetic field (in strength and/or 
direction) might redistribute the synchrotron emission power 
from the larger to the smaller angular scales, creating fake 
structures. The eff'ect is more relevant at lower frequencies, 
where the APS slope is expected to become flatter. We have 
verified this interpretation by performing a toy- simulation of 
pure Faraday depolarization in a simple slab model. We note 
that in the better sampled patches the polarization vectors are 
well ordered and similarly distributed at 820 MHz and 1411 
MHz, which is consistent with the slab model assumptions at 
least for the considered angular scales. We reproduced the ef- 
fects of Faraday depolarization on simulated polarized intensity 
maps, representing the intrinsic Galactic synchrotron emission 
at 408 MHz and 1420 MHz. We constructed an RM map inter- 
polating the data of Spoelstra (1984) and applied it to transform 
the input maps according to Faraday depolarization formulae. 
At 1420 MHz the input and output maps are almost identi- 
cal, whereas at 408 MHz the large-scale structure of the out- 
put map appears depressed and the small scale structures more 



abundant. Consequently, the corresponding 408 MHz APSs are 
flatter than the original ones (see Sect. 7.2 for a more precise 
description of the simulation process). It is therefore prefer- 
able to extrapolate the results (spectral index and amplitude of 
the APS) obtained at 1411 MHz (significantly less aff'ected by 
Faraday depolarization eff'ect than those at 408 MHz) to esti- 
mate the synchrotron APS at microwave frequencies. 

By interpreting the ratios of the APS amplitudes at 820 and 
1411 MHz (the latter smoothed to the lower angular resolu- 
tion of the former) in terms of synchrotron emission depolar- 
ized by Faraday eff'ect, we have identified possible RM ranges 
(see Sect. 7.1). This seems a reasonable approach at least in 
the case of patch 1, 2, and 3, since the polarization vectors are 
mostly aligned, indicating rather homogeneous physical condi- 
tions. Furthermore, taking into account the upper limit to RM 
derived from the available information on the degree of polar- 
ization at 1.4 GHz and considering that the maximum theo- 
retical degree of synchrotron polarization is ~ 75%, we could 
break or, at least, reduce the degeneracy between the identified 
RM intervals. For the three better sampled patches, the most 
reasonable value of RM are ~ 9 - 17 rad/m^. However, given 
the uncertainty on the measured polarization degree, RM val- 
ues in the interval - 53-59 rad/m^ cannot be excluded. Higher 
frequencies observations (at v ^ 2-3 GHz) would be extremely 
useful to break the degeneracy left by the available ~ 1 GHz 
data. 
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